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1. Introduction 

As a non-linear theory of pions, the Skyrme model |I| provides an approximate descrip- 
tion of hadronic physics in the low-energy limit. In this theory, the nucleon emerges 
as a non-perturbative solution of the field equations, or more precisely as a topological 
soliton. This model is also seen as a prototype which might be applicable in various 
physical contexts where one could expect soliton solutions to occur (e.g. condensed 
matter (baby Skyrmions), wrapped branes ...). More recently, this picture has re- 
gained attention since it could provide an explanation for the newly discovered hadronic 
states H, i. 

The original Skyrme lagrangian is a naive extension of the non-linear sigma model 
consisting of a fourth-order field derivatives term. This is nonetheless sufficient to 
stabilize the soliton against scale transformations and to reach at least a 30% accuracy 
with respect to physical observables. In order to incorporate effects due to higher-spin 
mesons and improve the fit on most observables a number of alternative Skyrme-like 
models which preserved the form of the original lagrangian while extending it to higher 
orders has been proposed and analyzed 0, |, 0, H • Unfortunately, all those models are 
limited in that they do not admit exact analytic solutions. Indeed, very few analytic 
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soliton solutions are known namely the one- dimensional sine-Gordon equation , the 
KdV equation, the instanton and some other special cases. This is because a necessary 
condition for a soliton-like solution to exist makes the search for an analytic solution 
rather complex and its discovery accidental. 

In the absence of exact analytical solutions, the only alternative to numerical treat- 
ment is the use of aptly chosen analytical forms which provide sometimes a reason- 
able approximation but which may not reproduce the correct behavior in the limits 
r 0, oo. Apart from greatly simplifying calculations of physical quantities, a great 
deal of information can be extracted from such analytic forms. For instance, symmetries 
and general behavior of the solution are much easier to analyze, and the characteristics 
of each Skyrme-like model also become more explicit. An analytic form would also 
prove useful in the stability analysis of the soliton, both classical and quantum, and 
in the calculation of multi-Skyrmion interactions. For example, one can analyze the 
quantum behavior of the Skyrme model soliton based on a family of trial functions, 
taking into account breathing motions and spin-isospin rotations or use these solutions 
to examine the two Skyrmion interactions. For a more accurate analysis, one has to 
resort to numerical computations which can be time consuming (e.g. the full numerical 
solutions for lower topological (or baryon number) < 1 7 of Houghton et al ) . 

The purpose of this work is to find multi-skyrmion solutions analytically assuming 
rational maps. For the sake of simplicity, we consider a class of models which is at 



most of order six in derivative of the fields. Floratos and Piette |jTO[ have indeed shown 
that for A^ < 5, the solutions have the same symmetries as the pure Skyrme model 
which are well represented by the rational map ansatz. The general character of the 
calculation also allows to extract an exact solution for the A^ = 1 hedgehog solution 
or for convenience in some cases, simplified approximations of the solutions for any 
A^. We proceed in three steps: First, we write the series expansion of the solution 
near r = and find the recursion relation for the coefficient of the series (Section III). 
Second, an analogous analysis is performed in the limit r — > oo. Finally, we propose a 
parametrization for the full solution and use the former recursion relations to set the 
parameters (Section IV). The calculations are somewhat intricate but once the relations 
are found they can be easily implemented as a computer algorithm both in numerical 
or in symbolic calculations. The method described here is also versatile in the sense 
that one could propose a different or more appropriate parametrization and still use the 
same recursion relations for r — ^ 0, oo to fix the parameters. We shall conclude with 
comments on the advantages and limitations of the procedure and on various ways of 
improving the convergence of the series (Section VI) where the chiral angle shows an 
abrupt behavior (e.g. larger A^). 
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2. Skyrme model 



The Skyrme model is defined by the lagrangian: 

1 

C = ^ Tr L.L^ + Tr f^.J'"' (2.1) 

16 ^ 32e2 •'^ ^ ^ 

where = Wd^U and /^j, = with U = U{x) is an element of the SU{2) 

group. The first term, £i, corresponds to the lagrangian of the non-linear a-model. 
This becomes obvious when one substitutes the degrees of freedom in U by a- and 
TT-fields using U = -^{a + ir ■ n) where r are the three Pauli matrices. The second 
term, C2, contains four derivatives of the pion field and can account for nucleon-nucleon 
interactions via pion exchange. It was first introduced by Skyrme in order to prevent 
solitonic solutions arising in the non-linear a-model from shrinking to zero and thus 
allowing for stable topological solitons. 

The coefficients and e are respectively the pion decay constant (186 MeV) and 
a dimensionless coupling often called the Skyrme parameter. Here, we shall use more 
appropriate units and rescale the lagrangian according to 

(-^TrL.L'^) +i(lTr/,.r^) (2.2) 

With this normalization, lengths can be understood as units of while energy or 
mass as units of 

The lowest energy soliton was found by Skyrme himself and it takes the form 

[/(r) =exp[iT-?F(r)] (2.3) 

where F{r) is called the chiral angle or profile function of the solution, and r is radial 
unit vector. Its spherical symmetry and the hairy configuration of the spin and isospin 
pointing out at infinity has earned it the name of hedgehog solution but more techni- 
cally, such field configuration constitutes a map from physical space onto the group 
manifold SU{2) and is assumed to go to the trivial vacuum for asymptotically large 
distances. The latter constraint allows imposing that U{r —>■ 00) — > 1 from which one 
may derive the existence of a topological invariant associated with the mapping. The 
originality of Skyrme's idea was to identify this invariant, i.e. the winding number N, 

N = (factor) j d^rBo with = e^^p^L^LfL" 

where Bq is the topological charge density, with the baryon number. 
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For the static hedgehog configuration ( p.3|) , the energy density is the sum of 



Si = -^Ti LiV = [2a + b] (2.4) 



£2 = -—Trfi,r = a[a + 2b]. (2.5) 



1 

16 

with a = ^^^-^ and b = F'"^. Notice that although the second term in ( ^.2| ) is quartic 
in the derivatives of the pion field, the corresponding energy 82 remains only quadratic 
in F'. 

Despite its relative simplicity, the Skyrme model has deep implications and is 
rather successful in describing low-energy hadron physics. Yet, it cannot be consid- 
ered seriously as potential candidate for the full low-energy effective theory of QCD. 
For example, there is no compelling reason (other than simplicity) and certainly, no 
physical grounds to exclude higher-order derivatives in the pion field from the effective 
lagrangian. On the contrary, large analysis suggests that the bosonization of QCD 
should involve an infinite number of mesons, which implies that in the decoupling limit 
(or large mass limit) for higher spin mesons that it leads to an all-orders lagrangian for 
pions. Several attempts to construct a more realistic effective lagrangian were made 
either by adding vector mesons to the Skyrme picture or higher-order derivatives terms 
to the lagrangian 

In the latter approach, one adds a sixth-order term involving the baryonic density 

B, i, 

Cj = CjTr [5^5^] = 3a% with = e^.p^L'^L^L'^ 

where cj is a constant. The term £3 = -^Ti f^^i^f^^f^ ^ proposed in @] leads to an 
identical contribution to the static energy density,. 

£3 = ^o?b (2.6) 

Allowing for all-orders in derivatives of the pion fields, one can demonstrate that models 
can be constructed in terms of £1, C2 and £3 alone and that in a rather simple class 
of these models the static energy density is a combination of the three invariants, 
81,82 and £^3, which remains quadratic in F' and thus tractable. Moreover the full 
lagrangian in this case can be easily written using a generating function [|ll[ . 



The purpose of the present work is to find multi-skyrmion solutions analytically 
and, for the sake of simplicity, we shall only consider a subclass of these models, those 
that can be written as a linear combination of £1, C2 and £3 i.e. the most general 
sixth-order lagrangian. The static energy density for the = 1 hedgehog ansatz in 
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such a generalized model is simply given by 

3 

£" = ^ hmSm = 3x(a) + (6 - a)x'(a) 

m=l 

where hm are coefficients and x{^) = Ylm=i ^mx'^- Integrating over volume leads to 
the mass of the soliton 

Ms = 4:TT r^dr[3xia) + {b- a)x'ia)] 
Jo 

Using the same notation, the chiral equation becomes: 



= x'ia) 



F' sin F cos F 
F" + 2— - 2 



+ ax (a) 



F' .n, COS F sin F cos F 

-2- + F'2— — + 

r sm r 



with a 



sin^ _F 



The Skyrme lagrangian corresponds to the case x{(^) 



Xs[a) 



N > 1 Skyrmions and rational maps 

For N > 1 Skyrmions, we shall assume that they are conveniently described by ap- 
propriate rational conformal transformations on the hedgehog solution as suggested 
by Houghton et al. 0. This assertion has been verified numerically up to iV = 5 by 
Floratos and Piette |10[ where they find that the solutions exhibit the same symmetries 



as the pure Skyrme model. 

The rational map ansatz lies on a naive connection between rational maps which are 
maps from S"^ h-^ S"^ and Skyrmions which are maps from i-^ 5^. It is then possible 
to identify the domain S*^ with concentric spheres in M^, and the target S*^ with spheres 
of latitude on S^. Amazingly, these are accurate to a few percent with respect to lowest 



energy solutions obtained through lengthy numerical calculations. loannidou et al. 
later generalized the ansatz in the context of SU{N) Skyrme models using harmonic 
maps from S"^ into CP^^^ . 

Rational maps are usually described in terms of the complex coordinate z = 
tan(^/2)e*'^ which correspond to stereographic projections on conventional polar co- 
ordinates 6 and ip. The point z identifies a unit vector 

1 



1 + \z\^ 

A rational map is a conformal map defined by 

q{z) 



2n 



5 



where one of the polynomials p{z) and q{z) is at least of degree N. One can then 
associate a unit vector to the rational map R{z) as follows: 



1 



1 + |i?P 



{R + R,R- R,l- 



In the context of the Skyrme model, the rational map ansatz propose a solution of 
the form 

f/(r, z) = exp(i ur ■ TF{r)) 

which reduces to the hedghog ansatz (|2.3|) when = 1. With boundary conditions 
F(0) = TT and F{oo) = 0, the baryon number turns out to be A^, the degree of R{z). 
Using this ansatz, the static energy density for this class of models takes the form 

oo 
m=l 

where now = F'"^ = b and = ap{z) with 



1 + 


z 


2 


dR 


1 + 


R 


2 


dz 



A remarkable advantage of the rational maps lies in the fact that the separation of the 
angular and radial dependence of the solution is preserved. Integrating over to obtain 
the mass of the soliton, we have 



M.s 



/ r'^dr / dVt {ap)^^^ [(3 — m)ap + mb\ 

J m=l 



where the integration over solid angle reads 



dVt 



2i dzdz 



[l + \z 



2\2 



and the factor 2i dzdz /{I + is equivalent to the usual area element on a 2-sphere 
sm6d6d(p. 
Defining 

2m 



~ ' Att^ ~ An 



1 + 


z 


2 


dR 


1 + 


R 


2 


dz 



2i dzdz 



1 + b 



2^2 ' 



and 



E 



m=l 



(y.rn.Q' 



X2 



(3.1) 



m=l 
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where am = hmlm^ = hmlZ-i^ write the expression for multi-Skyrmion 

masses 



' r^dr (3xi(a) - ax'i(a) + 6x2(«)) 





(3.2) 



Some of the angular integrations are trivial for rational maps: = 1, = N, = 1. 
Of course, the solutions depend on the model, i.e. the weight of each terms in xi and 
X2, or more precisely on the value of the each coefficient h^- This in turn determines 
which rational maps R, and values of angular integrations I^, minimizes the mass. 

Prom now on, we shall restrict our analysis to models of order six in derivatives of 
the pion field by setting coefficients /i^ = for m > 3. These models correspond to the 
most general case where the choice of rational maps affects a single non-trivial angular 
integration and the expression for the mass simplifies according to 



Xi{a) = hiNa + /la/s^a^ 

X2{a) = hia + h2Na^ + h^I^a^. 

Using the Derrick arguments, a stable soliton exist for = 1 only if is positive. This 
term prevents the Skyrmion from shrinking to zero size against a scale transformation. 
In this case, to find the minimal energy configuration at fixed (fixed degree for R{z)) 
one proceeds as follows: (i) minimize as a function of the coefficients of polynomials 
p{z) and q{z) and (ii) find the profile function F{r) which minimizes the energy. Note 
that since is positive and can go up to infinity ioT N > 1, negative values /i2 are 
physically excluded, as minimization would lead to an infinitely large negative soliton 
mass. One then concludes that rational map configurations and symmetries which 
minimizes the mass for a general sixth-order model (positive hi,h2 and h^) are the 
same as for the Skyrme model since one must simply minimize in both cases. 
The chiral equation becomes accordingly: 



O^hi (f" + 2y]+2N[ -hi 



+ L 



N 



-2ho 



sin^F 



sin F cos F 



sin F cos F 



+ 3/i3 



+ h2 



sin^F 



sin^F 



F" + F 



1/2 



F'' 

F" + 2— 
r 



cosF 
sinF 



(3.3) 



Our goal is to obtain an analytic expression for the solution of the chiral equation. 
For this purpose, we write the chiral equation in a form that is more convenient using 
(j) — cos F : 
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= -2rVi(l-0')V + Xi(l-0')'0 (3.4) 

+ X2 (-0" (1 - 0') - i<P'? <pr' - 2r (1 - 0^) 0') + X2 {^r (l - 0^) 0' + (0')') (l - 0^) 



(3.5) 



where Xi = Xi 



We would hke to represent its solution 0(r) in the form of a series. Its explicit form 
is a matter of choice as long as it reproduces the solution adequately over all values of 
r which requires the convergence of the series. Otherwise, what choice of series is best 
would have to be judged on a number of criteria which have to be met at least minimally: 
(1) rapid convergence, i.e. fewer terms are necessary to reach numerical precision with 
respect to the exact solution, (2) mathematical tractability e.g. recursion formula for 
the coefficients of the series can be written and (3) other aesthetical criteria, e.g. a 
form of series leading to an energy density and other physically quantity that can be 
integrated analytically. 

Solutions can be constructed simply by adding terms to approximate solutions 



proposed in the past |T3|, [Tj, [T^, [18], |T9|. However, a quick analysis reveals 

that some of these candidates must be discarded since they do not behave properly at 
small or large distances. In this work, we propose a solution which implements these 
constraints and a procedure to allow for the convergence of the series. However, the 
first step of our analysis will be independent of a specific choice of series representation 
for the solution. It consists in finding the small (r 0) and large (r — >■ 00) distance 
behavior of the solution as a series expansion in powers r and 1/r respectively. This is 
computed using the chiral equation. Once the small and large r expansion are known 
explicitly, the information will serve to propose an appropriate form for the series and 
ultimately to calculate the coefficients of the series. 



4. Small distance behavior 

First, let us write as a power expansion in r in the limit r 0. The chiral equation 
( |3.5|) is symmetric with respect to the change r —>■ — r, 0(r) 0(— ?") = 4>{r), which 
simplifies the power expansion since only even powers of r contribute: 

00 

0(r) = ^a„r2'^ (4.1) 

n=0 

with boundary condition ag = —1 since 0(0) = —1. 
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For convenience, let us write explicitly the functions xi and X2 in the chiral equation 

m 



oo 



= r V XI - 3)a„a"+^ + (0')^ J] n{n - 2)/3„a"-^ (4.2) 

n=0 n=0 

oo oo 

+ 2r0' J2 ^(^ - 2)/3„a" - r V" 

n=0 n=0 

Substituting by the series ( |4.1D and 0', 0", {(p' f and a" by the appropriate expres- 
sions 

in the last equation, we get 

3 oo oo 
n=l k=0 1=0 

The coefficients and c„_fc are given respectively by 

hi = iaoal 

bk = bk + Abk for A; > 1 
co,o = 1 

Co,fc = for > 

Cn,o = (— 2aoai)"' for n > 
Cn,fc = Cn,fc + Acn,fc for n. A; > 

with 

k—l k—l /k—m 

bk = 4ao X + 1 - q)agak+i-q + 4 ^ I ^ g(A; + 1 - m - g)agafc+i_,„-g 

g=2 m=l \ q=l 

/k—l k q 

Cn,k = ^ ( ^{q{n + 1) - k)2aoaq+iCn,k-q + ^{q{n + 1) - A;) ^ a^ag+i.^c^^fc-g 
Ac„,fc = -2"naoa"~^afc+i 



g=l g=l r=l 



where for later convenience, we have isolated the contribution proportional to the co- 
efficients Ofc with highest k in the terms Afo^ and Ac„,fc respectively. 
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The equation holds for arbitrary r, which imphes that it holds for any arbitrary 
power of r. Isolating the term in r^*^*^"^-' we get an expression of the form 

= Ak-iak + Bk-i. (4.3) 

where Ak-i and Bk-i are respectively 

3 

Ak-i = XI 2^ {n - 3) (n + 1) + /5„ {n^ - 2n - 2 - 2^ + k)) , 

n=l 



k-1 



n=l m=2 



m 

+ 2m {2n — 3 — 2m) amCn,fc-m + 2(2n — 5)aiC„ ^^i). 

From inspection we see that both Ak-i and -Bfc-i depend on lower- index coefficients am 
(i.e. m < fc) so we can write (|4.3| ) in the form of a recursion formula for the coefficients 
ttk for k > 1, 

ttfc = — T~~~" with ttQ = —1. 

The formula requires a single input parameter ai (which must be positive), aside 
from the boundary condition oq = —1. cti is related to the slope of the profile function 
at r = as follows 

F{r^0) = ^T-^/2a^r + O{r^) or 0(r ^ 0) = -1 + Oir^ + O(r^). 

It is also easy to verify that ak depend on all a„'s and (3nS whose values are 
prescribed by both the model and the topological number of the solution (see (|3.1|) ). 
This has an interesting implication: The small r behavior of the solution is characterized 
by all the terms in the lagrangian and not only the one with highest-order derivative 
in the pion field as one might have suspected. 



5. Large distance behavior 

Now let us examine as a power expansion in 1/r in the limit r — oo Again since the 
chiral equation (|3.5|) is symmetric with respect to the change r — — r, 0(r) 4>{—r) = 
(f){r), the power expansion only contains even powers of 1/r. Accordingly, we write 

oo 

= 5^ a^r-^™ (5.1) 

m=0 
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Allowing for the boundary condition 0(oo) = 1, we set oq = 1. It is also easy to verify 
that 

di = and 0,3 = 0. 
According to (|5.1| ), we write the following expression as power series 



where the coefficients bk and Cn,k are given respectively by 

65 = 16000-2 

bk = bk + Abk for > 5 

Cn,0 = 

Cn,k = for < 3n 

Cn,3n ( ) 

Cn,k = Cn,k + Ac„_fc for k > 3n 

with 

^ k-4 k-1 k-3-m 

6fc = 4ao^/(A; - 1 - /)d/afc_i_; + 4^a„ ^ l{k - I - m - l)aiak-i-m-i 

1=3 m=2 1=2 

Abk = 16(/c - 3)aoa2dfc-3 

k—3n ^ fc— 3n— 1 

'^n,k = -n (-2aod2)""^ agak-3n+2-q - . 5-^ + 1) - + 3ra)ci,g+3C„,fc_g 

Acn,fc = -2ao^ (-2aoa2)"~^ dfc-3n+2- 

Here the contributions proportional to the coefficients with highest k are written 
exphcitly in Abk and Ac^^^ respectively. 

Substituting (p by the series (|0| ) 0', </>", </> (0')^ and a by the above expressions we 
arrive at 

00 00 fc 



= ^ ^ ^ ^n{n - 3)anamCn+i,k-mr ^''"^^ + n{n - 2)PJ)mCn-i,k~mr "^^ (5.2) 

(5.3) 



n=0 A;=0 m=0 

00 

2k 



-2n{n - 2)Pn2mdmCn,k-mr - ^ ra/5„2m(2m + l)a„c„,fc-m?" 

n=0 

Again, since the equation must hold for arbitrary r, it implies that it holds for an 
arbitrary power of r. Therefore isolating the term in r~^''~^^ leads to the relation 

= Ak-iak + Bk-i (5.4) 
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with 



Ak-i = 4a2 (-4ai + /?i {2k^ - 5k + 6)) 



k-2 

OmC2,A; — m+4 j 



m=l 

fc+2 k-l 



\ m=0 m=3 

oo fc+3 

+ ^ ^ (^?^(?^ - 3)a„amCn+i,fc-m+4 + n{n - 2)/3„6,„c„_i,fc_m+3 

n=2 m=0 



-2n (2n - 3 + 2m) m/3„amC„,A:_m+3) . 

After inspection, we find that the right hand side of equation (|5.4|) is trivially zero for 
A; < 3 and but otherwise contains only terms in with m < k. Furthermore, Ak-i 
and only depend on lower-index coefficients {m < k — 1). This allows writing 
a recursion formula for the coefficients for A; > 3, 

Bk-i 

ak -- 



A 



k-l 



Recalling that Q,Q — 1, 0-1 = and as = 0, the only remaining unknown parameter is 0,2 
which must be negative. This latter parameter fixes the dominant contribution to the 
profile function in the large r limit. 



F(r^^) = ^ + 0(i) or ^ ^) = 1 + | + + ... 

Moreover, the first coefficients in the large r expansion depend only on the lowest-n 
a„'s and /5„'s. This means that the large-r behavior is not very sensitive to the model 
and to the topological sector of the solution. In other words, higher-order derivative 
terms and the topological sector begin to contribute only when the sub-dominant terms 
in r~^*^ become important. 



6. Solutions for Skymions 

We have found recursion relations for the coefficients of a series expansion in powers 
r and 1/r for small (r 0) and large (r —>■ oo) distances respectively. We may now 
proceed to construct a full solution. The explicit form of the series is to a certain 
extent a matter of choice. Of course, it should represent the solution adequately which 
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means that convergence is a prerequisite and that small (r — 0) and large (r oo) 
distance behaviors must be reproduced. Other criteria may also be suggested such as: 
(1) faster convergence, i.e. fewer terms are necessary to reach numerical precision with 
respect to the exact solution, (2) mathematical tractability e.g. recursion formula for 
the coefficients of the series can be written, (3) the energy density can be integrated 
analytically, etc... 

We propose here a solution and a procedure similar to that in but which applies 
to 4> = cosF instead of F. Taking in account the small (r — > 0) and large (r —>■ oo) 
distance behavior, we can write recursion relations for the coefficients of the series. We 
shall see that this particular ansatz also allows in principle, to calculate the mass of 
the soliton analytically which was not possible in [T^. The solution has the form 



Mr) = = 1 + E 2 2%m+3 (C2™(r' + rl) + c^m+ivl) (6.1) 

^0 m=0 + ^o) 

The small distance expansion is given by 



+ {2q - k + 4:){2q - k + 3) (c2g-2fc+2 + C2g-2fc+3)) 

^ J + (3co + 6ci - 4C2 - 5C3 + C4 + Cs) 

Matching the coefficients at small distance with those of (^4.3|) , we get 

ao = 1 + Co + ci = -1 

and for k > 1 

, ^.k {k + 2){k + l){-l)\ {2k-r)\{-iy 

{Tq) Qk = (ao - 1) + C2fc-2 + C2k + C2fc+1 + ^ 



2 v-^u . , , ^z«+. , ^H(2fc-2r + 2)! 

■ {{2k - 2r + 2)(2fc - 2r + l)c2fc_2r-2 + (2A; - r + 2)(2/c - r + 1) (c2fc^2r + C2fc-2r 
We invert the relation to get the odd-index coefficients and obtain 

ci = ao — 1 — Co = — 2 — cq (6.2) 
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and for /c > 1 



C2k+i = [ro) ak-C2k-C2k-2 7^ (ao - 1) (6.3) 



k-l 



r=l ^ ' 

+ {2k-r + 2){2k - r + 1) {C2k-2r + C2k-2r+l))- (6.5) 

On the other hand, the large distance behavior reads 

-2 oo / 2\ -9-3 



q=0 



+1> 



(2g- A; + 2)! (-!)''■ , , , , ^ 
+ g (fc + l)!(2,- 2fc + 2)! 1) + 

Again by inspection 

— 2co + ci + C2 = or C2 = 2co — ci = 3co + 2 

and matching the coefficients at large distance, we get the even-index coefficients of the 
solution which then reads 

/ 2\-fc-2 (2A; - r)! (-1)'' , , , , ^^ ^ 

C2k = ak+2 (ro) - 2^ (r + l)\{2k - 2r)\ ^ ^ ^ ''^ C2fc_2r-2 + (r + 1) C2k-2r-l) ■ 

(6.6) 

Inserting (|67^), ( |6^ ) and ( |6.6| ) in (|6.1|), we end up with a solution which is now 
completely determined by three parameters, ai, 02 and tq. The ffist two parameters 
tti, 0,2 depend on the behavior near r ^ and r —>■ 00 respectively whereas tq can be 
interpreted as an intermediate scale. The strategy adopted here to reach the solution 
lies upon the fact that for a given k, the even-index coefficients C2k in (|6.6|) depend 
on lower-index coefficients with m < 2k and 0^+2 which is fixed by the {k + 2)*^ 
coefficient of the large r expansion or ultimately 0,2. On the other hand, the odd-index 
coefficients C2k+i in ( |6.4|) requires Cm with m < 2A; + 1 and au which is determined by 
the behavior of solution at small r or the parameter ai. This has two effects: First, the 
procedure refines the solution by alternatively the matching of the odd- and even-index 
coefficients, which requires higher and higher derivatives of the solution at r — > and 
r —>■ 00 respectively ( |6.1D . Secondly, since it relies on matching the coefficients and 
dk+2, the accuracy of the series and convergence must improve near the end points as 
k increases. 
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It remains that the parameters ai, 0,2 and tq are still unknown at this point of 
the procedure. Several approaches to find these parameters are possible. One such 
procedure consists in introducing three more constraints on 0^^. For example, we can 
use the chiral equation ( p.5|) at three intermediate points to set the values of ai, 02 and 
tq. We could also choose to impose a continuity condition at a given r = tq assuming 
that the series ( |4.1| ) and ( |5.1| ) apply to r < tq and r > vq respectively. This procedure 
would lead to a solution which has the virtue of being completely determined by the 
equation of motion, but they do not guarantee that the set of parameters would be 
the best to render the exact solution. Since these methods would require numerical 
calculations at some point anyhow, we adopt a more practical approach which consists 
of finding numerically the set of values for ai, 0,2 and ro which minimize the mass of 
the soliton. 

In principle, the series ( |6.1|) contains an infinite number of terms which allow in the 
end to reach the exact the solution. For computational reason the series is truncated 
which means that, the number of terms in the series will affect not only the precision 
of the series but also the values of the parameters Oi, 02 and rg that minimize the mass 
of the soliton. The larger number of terms in the series, the closer we get to the exact 
solution and values of the parameters ai, 02 and tq. 

We could attempt to compute the mass of the soliton analytically. Starting from 
(|3.2|) we have 

Ms = 4:TT / 2u^ (1 - u)^ du {hi {2Na + b) + {l^a + 2Nb) + Shl^a^b) 
Jo 

where we changed variable r to r ^ -u = ^2^^^^ or r = rQ^^J^. Substituting (p by our 
solution 

K 

= 1 + J] n*^(l-^z)'+'(c2fc + C2fc+i(l-^i) ) (6.7) 

k=0 

where K is finite but sufficiently large for accuracy, and assuming that we can expand 
the integrand in powers u and (1 — m) , it becomes possible to integrate the expression 
analytically since all integrals can be cast in the form of 

Mab = C W'-Hl - uf-^du = ^}r^^ for Re(a), Re(6) > 0. 
Jo r(a + b) 

Unfortunately, even for the simplest case, the Skyrme model, the calculation becomes 
inefficient and quite impractical as K increases. So we resort to numerical integration 
which in this case is much faster and proceed to minimize the mass in terms of the 



15 



parameters ai, a2 and tq . As an example, we show the result of the first few values of 
K for the Skyrme model,, for = 1, which requires hi = h2, = and ll = I2 = 1 





K = 1 


K = 2 


K = 3 


K = A 


Num. 


Mass 


1.23268 


1.23174 


1.23151 


1.23151 


1.23145 


ai 


2.17741 


2.05754 


2.02667 


1.99411 


2.01508 


0,2 


-2.39428 


-1.62650 


-1.22111 


-3.15231 


-2.33204 


ro 


1.05398 


1.12095 


1.18501 


0.939329 





The exact mass is obtained within an accuracy of 0.005% for K as low as 3. Note 
however that the convergence of ai and 02 towards their exact values is not as efficient. 

7. Advantages and limitations of the approach 

An extensive comparison of the approximate solution ( |6.7D with respect to the exact 
numerical solution suggests that the reliability of the approach depends largely on the 
model considered and topological sector. The recursion formulae for the Cfc coefficients 
ensure a perfect agreement in both limits r — and r — > 00. The agreement is still 
preserved when the solution remains smooth for all values of r, as in the case of low 
A^ soliton in the Skyrme Model. Nevertheless, discrepancies between analytical and 
numerical results begin to appear when we consider sixth-order models or when the 
topological number A^ increases. The solutions in these cases are characterized by a 
sharp behavior of F near F{r) = 7t/2 (or near 0(r) = 0). 

In order to understand the origin of these discrepancies, it is instructive to look 
at the quantity a = sin^ F for small, large and intermediate r. For the interme- 
diate region where sin^ F reaches its maximum, the energy density is found to be the 
largest and F is almost linear. On the other hand, F" and of course cosF are rela- 
tively small. Therefore, which term dominates the chiral equation (|3.3| ) in that region 
depends on the relative weight of the coefficients hi, 2Nhi,2Nh2, /^/i2 and /^/is and 
on the highest power of a. As A^ increases, increases approximately as A^^ and the 
term proportional to /is/^ is expected to dominate if /12 and hs are of the order same of 
magnitude. This implies that most of the variation of F should occur in the interme- 
diate region since the dominant term is proportional to which is suppressed outside 
this region. This leads to a configuration of the energy density which is localized on 
a shell of decreasing thickness as A^ increases. Unfortunately, the recursion formulae 
for the Cfc coefficients only apply for the end regions r and r —>■ 00 which do not 
contribute much to the soliton mass. So we can expect the approximate solution ( |6.7|) 
to lose its accuracy when A^ gets larger. In other word, although accuracy is expected 
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to increase in the r — > and r — > oo region, it does not improve in the intermediate 
region where most of the energy is concentrated. 

Various ways of improving the analytical solution may be considered. For instance, 
several other trial functions have been used in conjunction with the small and large 
r recursion formulae, including polynomials series and Pade approximants, but we 
found no general form that would eliminate the aforementioned discrepancies. In fact, 
there is no reason a priori why any trial function should evade the conclusion of the 
last paragraph. A better improvement would fix the behavior of the profile angle and 
describe accurately the energy density in the intermediate region where its contribution 
is more important for any N. This could proceed through the construction of a recursion 
formula at an intermediate point and eventually its implementation into a series but, 
it implies rather complex intricate calculations and has not been attempted yet. 

Even though imperfect, we stress again that this approach is still probably the best 
alternative to complete numerical treatment, and that it can prove very useful whenever 
an analytical form of the solution is required, as we have already discussed regarding 
stability analysis. Because of the general character of the calculation, it easy find the 
exact solution for the N — 1 hedgehog solution for the pure Skyrme model. It is also 
possible — if one chooses not use the full solution — to extract approximate solutions 
for any since the first few terms of the series represent fairly good approximations. 
Moreover, this work also provide the first construction of the recursion formulae for 
the series expansion in the r — > and r — > oo limits which may be helpful in many 
calculations. 

This research was supported by the Natural Science and Engineering Research 
Council of Canada. 
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